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ABSTRACT 

We report on long-term monitoring of three anomalous X-ray pulsars using 
the Rossi X-ray Timing Explorer (RXTE). We present a phase-coherent tim- 
ing ephemeris for 4U 0142+61, and show that it has rotated with high stability 
over 4.4 yr, with RMS phase deviations of 7% of the pulse period from a sim- 
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pie fit including only v and v. We report on the continued timing stability 

of IE 2259.1+586, for which phase coherence has now been maintained over 
4.5 yr, as well as on the detection of a significant i) in 1.4 yr of monitoring of 
RXS J170849. 0-400910, consistent with recovery following a glitch. We note a 
correlation in which timing stability in AXPs decreases with increasing z>. The 
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timing stability of soft gamma repeaters in quiescence is consistent with this 

^ . trend, given their large spin-down rates. This trend is similar to one seen in 

radio pulsars, suggesting a connection between the three populations. We find 

no large variability in pulse morphology as a function of time. We present high 

signal-to-noise ratio average pulse profiles for each AXP, and consider them as a 

function of energy. We find a variety of different behaviors, and consider possible 

trends in the data. We also find no large variations in pulsed flux, and set ler 

"£o . upper limits of ~20-30% of the mean. 

d • 

Subject headings: pulsars: general — pulsars: individual (4U 0142+61, 
RXS J170849.0-400910, IE 2259.1+586) — X-rays: general 

W. 

1. INTRODUCTION 

Anomalous X-ray pulsars are an unusual class of astrophysical objects. There are cur- 
rently only five confirmed AXPs: 4U 0142+61, IE 1048.1-5937, IE 1841-045, RXS J170849.0- 
400910 (hereafter 1RXS 1708-4009), and IE 2259.1+586. All five are located in the Galactic 
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plane and two (IE 2259.1+586 and IE 1841-045) are coincident with supernova remnants. 
A sixth AXP candidate, AX J1845. 0-0258, is also coincident with a supernova remnant 
(Gaensler et al. 1999). 

The observational properties of AXPs (see Israel et al. 2002, for a review) are gener- 
ally quite different from those of conventional accreting binary X-ray pulsars (see Bildsten 
et al. 1997, and references therein for a compendium of accreting X-ray pulsar properties). 
In particular, AXPs have spin periods in the narrow 6-12 s range, while accreting X-ray 
pulsars have periods that span 0.002-1400 s. AXPs have luminosities in the narrow range 
10 34 — 10 36 erg s _1 . By comparison, typical accreting X-ray pulsars are extremely variable, 
with luminosities spanning a much broader range (10 34 — 10 39 erg s _1 ). AXP spectra are con- 
siderably softer than those of the typical accreting pulsar. Accreting pulsar spectra are gen- 
erally well described as a power law with a cut-off in the 10-20 keV range. AXP spectra, by 
contrast, are generally characterized by a blackbody spectrum of energy kT ~ 0.35 — 0.6 keV, 
with a hard energy excess that can be characterized by a steep power law of index —2 to —4. 
The timing properties of AXPs also contrast sharply with those of most accreting pulsars. 
All of the known AXPs, as we show in this paper, appear to be undergoing steady, prolonged 
spin-down, though some are more steady than others. By comparison, most accreting pul- 
sars are spinning up (on average), or alternating between spin- up and spin-down, and have 
poor rotational stability. 

Most importantly, the AXPs show no evidence for a binary companion. Specifically, the 
limits on the X-ray/optical flux ratio in AXPs rule out the presence of massive companions as 
in the conventional high-mass X-ray binaries (e.g. Mereghetti et al. 1992; Corbet & Mihara 
1997; Israel et al. 1999a). In addition, careful pulse timing of AXPs has failed to detect any 
evidence for binary motion of the neutron star on time scales of a few minutes to several days 
(Mereghetti et al. 1998; Wilson et al. 1999). The limits imply that only extremely low-mass 
companions are allowed in these systems, unless we are observing all AXPs face-on, which 
is unlikely. Furthermore, mass-transferring low- mass companions are also unlikely because 
the supernova remnant associations imply youth, and because these systems rarely survive 
the supernova. 

Currently, models for AXPs fall into two distinct categories. One class of models is 
that AXPs are young, isolated, ultra-magnetized neutron stars or "magnetars" (Thompson 
& Duncan 1996; Heyl & Hernquist 1997). If the spin-down of the pulsar is primarily due to 
magnetic dipole radiation, then the AXPs have enormous surface magnetic dipolar fields, in 
the range B ~ 10 14 — 10 15 G. The identification with magnetars is further strongly motivated 
by the similarity of the AXP emission to that of the soft gamma-ray repeaters (SGRs) in 
quiescence. Specifically, the latter have similar pulse periods, are spinning down (Kouveliotou 
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et al. 1998; Kouveliotou et al. 1999), and have X-ray spectra that are comparable to, though 
somewhat harder than, those of the AXPs, at least when not in outburst (Mereghetti et al. 
2000; Kulkarni et al. 2001). Independent evidence for the ultra-high magnetic fields exists 
in SGRs; for example, a ~ 10 15 G magnetic field is required to contain the radiation that is 
seen following major outbursts (Thompson & Duncan 1995; Thompson & Duncan 1996). 

The second category of model proposed to explain AXP emission is that they are accret- 
ing from a disk of material leftover from the supernova explosion that created the neutron 
star (Chatterjee et al. 2000; Alpar 1999; Marsden et al. 2001). Fall-back disks around neu- 
tron stars have been proposed in other contexts (Michel 1988); here, the AXPs represent a 
single, short-lived phase in the evolution of a neutron star having magnetic field only slightly 
larger than those of radio pulsars. In this phase, the neutron star is in a propeller mode in 
which the inner edge of the accretion radius lies between the magnetospheric and corotation 
radii. This model, in its original version (Perna et al. 2000), significantly overpredicted the 
expected optical/IR flux from the fall-back disk (Hulleman et al. 2000a,b; Kaplan et al. 
2001). A more recent version includes the effects of a thermal disk instability which may 
prevent the disk from expanding radially, hence allowing it to remain sufficiently dim to be 
consistent with the optical/IR limits (Menou et al. 2001). Nevertheless, difficulties remain 
with the accretion model. In particular, although the noise properties of a fall-back disk are 
currently unknown, the extremely steady spin down of two AXPs (Kaspi, Chakrabarty & 
Steinberger 1999) is at odds with that seen in most X-ray pulsars that definitely have ac- 
cretion disks (but see Baykal et al. 2001). Furthermore, very noisy timing behavior seen in 
one AXP (IE 1048.1-5937) shows no evidence for correlated flux variations, as are expected 
in fall-back disk accretion models (Kaspi et al. 2001). 

In this paper, we present a variety of observational results from our program of regular 
monitoring observations of AXPs with the Rossi X-ray Timing Explorer. We make use 
of dedicated as well as archival monitoring observations spanning several years for three 
AXPs. Initial results from this program were presented elsewhere (Kaspi et al. 1999, Kaspi, 
Lackey & Chakrabarty 2000) and focused on the timing properties of 1RXS 1708-4009 and 
IE 2259.1+586. Kaspi et al. (2001) and Gotthelf et al. (2001) present results for AXPs 
IE 1048.1-5937 and IE 1841-045, respectively. In this paper, we report for the first time on 
4U 0142+61, and update timing results for 1RXS 1708-4009 and IE 2259.1+586. In addition, 
we use the same data sets to study pulse morphology variability and energy dependence, 
as well as the stability of the pulsed flux. Our goal is to clarify important observational 
properties of AXPs in the hope of offering new, quantitative data for use in testing AXP 
models. 
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OBSERVATIONS 



The results presented here were obtained using the Proportional Counter Array (PCA; 
Jahoda et al. 1996) on board the Rossi X-ray Timing Explorer (RXTE). The PCA consists of 
an array of five collimated xenon/methane multi-anode proportional counter units operating 
in the 2-60 keV range, with a total effective area of approximately 6500 cm 2 and a field 
of view of ~ 1° FWHM. Our observations consist primarily of short snapshots taken on a 
monthly basis (see Table 1). In addition, we used a handful of archival observations; the 
exposures in these observations vary from very short (T ~ 1 ks) to very long (T ~ 20 — 40 ks; 
see Table 1). We used the GoodXenonwithPropane data mode, which records photon arrival 
times with l-/zs resolution and bins energies into one of 256 channels. Due to the soft 
spectrum of these sources, we analyzed only events from the top xenon layer of each PCU. 
Photon arrival times at each epoch were adjusted to the solar system barycenter, and binned 
with 31.25-ms time resolution. 



3. ANALYSIS AND RESULTS 

3.1. Phase Coherent Timing 

In the timing analysis, we included only those events having energies in a predetermined 
range so as to maximize the signal-to-noise ratio of the pulse. The energy cuts used were 2.5 - 
9.0 keV for 4U 0142+61 and IE 2259.1+586 and 2.2 - 5.5 keV for 1RXS 1708-4009. Each 
binned time series was epoch-folded using the best estimate frequency determined initially 
from either a periodogram or Fourier transform (though later folding was done using the 
timing ephemeris determined by maintaining phase coherence; see below). Resulting pulse 
profiles were cross-correlated in the Fourier domain with a high signal-to-noise template 
created by adding phase-aligned profiles from previous observations. The templates are 
shown in Figure 1 along with their respective Fourier transforms in Figure 2. The cross- 
correlation returns an average pulse time-of-arrival (TOA) for each observation corresponding 
to a fixed pulse phase. The pulse phase at any time t can be expressed as a Taylor 
expansion, 

0(t) = 0(t o ) + u (t - t ) + -u (t - t ) 2 + ^o(t - to) 3 + • • • , (1) 

where v = 1/P is the pulse frequency, v = du/dt, etc., and subscript '0' denotes a parameter 
evaluated at the reference epoch t = to- The TOAs were fit to the above polynomial using 
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the pulsar timing software package TEMPO 4 . Unambiguous pulse numbering is made possible 
by obtaining monitoring observations spaced so that the best-fit model parameters have a 
small enough uncertainty to allow prediction of the phase of the next observation to within 
~ 0.2. Typically this requires two closely spaced observations (within a few of hours of each 
other) followed by one spaced a few days later, and regular monitoring thereafter, as long as 
phase coherence can be maintained. 

To minimize use of telescope time while maximizing precision, our monitoring data 
consisted of frequent brief snapshots of each pulsar (Table 1). These snapshots suffice to 
measure TOAs to good precision. However, at each epoch, the period as determined by 
a Fourier Transform or epoch-folding has typical uncertainty of a few milliseconds. Thus, 
snapshots can determine spin parameters with high precision only when phase coherence can 
be maintained. 

For 4U 0142+61, we report here for the first time a phase-coherent timing solution that 
indicates that this AXP has been an extremely stable rotator over 4.4 yr of RXTE monitoring. 
However, a 2-yr gap in the data precludes unambiguous absolute pulse numbering; we find 
two equally viable phase-connections that differ by one pulse. Unambiguous phase count 
can be maintained in 1.3 and 1.4-yr segments at the beginning and end of the 4.4-yr span, 
respectively. Phase residuals after subtraction of spin-down models that include v and v only 
in the first and second segments have RMS of 2.7% and 1% of the pulse period, respectively, 
and are featureless. The best-fit model parameters for each segment, as well as the two 
indistinguishable models for the entire 4.4-yr span, are given in Table 2. Timing residuals 
for solution "A" (see Table 2) are shown in Figure 3; those for solution "B" are qualitatively 
similar. 

It is interesting to compare our results with frequency measurements made for this pulsar 
with various X-ray telescopes over the past ~20 yr. Figure 4 shows these data (taken from 
Israel et al. 1999b) along with our RXTE ephemeris. The top panel shows v versus time; 
the thick solid line is our phase-connected ephemeris (the distinction between A and B is not 
visible on this plot) over the 4.4-yr span RXTE has monitored it, while the thin line is the 
extrapolation of our ephemeris backward. The uncertainties are smaller than the thickness 
of the lines. The bottom panel shows the same data but with the linear and offset trend 
removed. All but one data point (taken with EXOSAT) agree with the back-extrapolation. 
That one point suggests that the pulsar may exhibit deviations from a simple spin-down law; 
these could be due to random noise processes or glitches. Further observations can test this. 

For 1RXS 1708-4009, phase coherent timing has been accomplished in the 1.4 yr since 
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the glitch reported by Kaspi et al. (2000). In these data, we find a significant positive v (see 
Table 3). This indicates a decay of the negative z>, expected for long-term glitch recovery, as 
seen in glitching radio pulsars (e.g. Shemar & Lyne 1996). Phase residuals after subtraction 
of the best fit model are shown in Figure 5. We note that the value for z> reported in Table 2 
is different by 3.8a from that reported by Kaspi et al. (2000). This is likely due to the 
unmodeled influence of v in the earlier, shorter data set. 

The rotational stability of IE 2259.1+586, first reported by Kaspi et al. (1999), has now 
persisted over 4.5 yr, although the inclusion of v has been necessary. Phase residuals after 
subtraction of a simple model that includes only u, v and i> have RMS under 1% of the pulse 
period (see Table 3, Fig. 5). We note the presence of very low- level systematic trends in the 
most recent data; their origin is unknown. 



3.2. Pulse Morphology Changes 

To search for pulse profile changes, the profiles were first phase aligned using the tem- 
plates and the same cross-correlation procedure used for timing. Each data profile (D) was 
fit to the template (T) by adjusting two parameters, T and DC, to minimize a y 2 statistic, 
where 

2 1 ^ (D t -(F-T t -VC)) 2 

Xv N,-A^ ^ + ^-4 • U 

Here N^ is the number phase bins, and Or> and o^ are the errors associated with D and T, 
respectively. The resulting data profile was subtracted from the template to yield "profile 
residuals" for each observation. 

To search for pulse profile variations, using the optimal JF and VC, we calculated the 
probability of each xl statistic given their expected distribution and the number of degrees 
of freedom [y = N^ — 4). The procedure was repeated for each pulsar with N^ = 8, 16, 32 in 
order to have sensitivity to a variety of types of pulse profile changes. We have not detected 
any large pulse profile variations. This justifies our other analysis procedures which assume 
a fixed profile. However, in a handful of observations, we have found evidence for low-level 
pulse profile changes, at the ~ 3a level for all three sources. These are intriguing; however 
at this level, they require confirmation using longer monitoring observations. 

It is difficult to set quantitative upper limits on the amplitude of pulse profile changes 
to which we were sensitive, as these depend on the shape of the change, and vary depending 
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on the length of the observation. Typically, RMS profile residuals are ~20% of the pulse 
peak, although this varied from 4-40%. 

Iwasawa et al. (1992) reported a significant change in the pulse morphology of IE 2259.1+586 
in 1.2-14 keV GINGA observations obtained in 1990, such that the leading pulse had am- 
plitude roughly half that of the trailing pulse (compare with Fig. 1C). If correct, this has 
important implications for the magnetar model, which only predicts such pulse morphology 
changes in the event of a restructuring of the magnetic field, as might occur following a 
major SGR-like outburst. Although it is difficult to quantify upper limits on the amplitude 
of pulse profile changes generally, we have done Monte Carlo simulations to see whether we 
were sensitive to the same change observed by Iwasawa et al. (1992). In particular, using the 
original 1990 GINGA pulse profile (kindly supplied to us by B. Paul and F. Nagase) with 
random Poisson-deviated noise added, we repeated our pulse change morphology analysis as 
described above. Our simulations show that we would have detected the change with high 
(>99%) confidence in 25% of our observations (those with integration times >10 ks), and 
with moderate (>90%) confidence in 45% of our observations (those with integration times 
of >7 ks; see Table 1). 



3.3. Pulse Profile Energy Dependence 

AXP spectra are generally best fit by a two-component model consisting of a photo- 
electrically absorbed blackbody with a hard power-law tail (Israel et al. 1999b). Whether 
these two components are physically distinct is an open question . To investigate this, we 
compared the pulse profile morphology in two energy bands, 2-4 keV and 6-8 keV. These 
were chosen because the 2-4 keV band has a significant blackbody component (~ 40% — 60% 
of total flux, depending on source and model) and the 6-8 keV band is greatly dominated 
by the power-law component (^ 95% of total flux), (see Ozel et al. 2001) 

Figure 7 display the average pulse profiles of 4U 0142+61, 1RXS 1708-4009 and IE 2259.1+586 
in the energy bands 2-4 keV and 6-8 keV. Note that the scaling in these plots was chosen 
to minimize the x 2 of the difference between the soft and hard profiles. Thus the only infor- 
mation that these plots convey is the relative amplitudes of the features of the profile. The 
minimization of the \ 2 statistic (Eq. 2) between profiles in the two energy bands was done 
in the same manner as for the pulse morphology changes (see §3.2). 

Another way of characterizing the pulse morphology as function of energy is to consider 
the pulse profile harmonic content in different energy bands. Figure 8 displays the power 
of the n harmonic (A n ) in units of total power in the remaining harmonics (A tota i = 



A\ + A2 + ■ ■ ■ + An^ , where N<f, is the number of phase bins), versus harmonic number n, for 
each source. 

From Figure 7 and 8 it is clear that all three of the AXP pulse profiles vary significantly 
as a function of energy, but to different degrees. 



3.4. Pulsed Flux Time series 

Given the large field-of-view of the PCA, the low count rates for the sources relative to 
the background, and the fact that, for example, IE 2259.1+586 is in a supernova remnant, 
total flux measurements are difficult with our RXTE data. Instead, we have determined the 
pulsed component of the flux, by using the off-pulse emission as a background estimator. 

Data from each observing epoch were folded at the expected pulse period as was done 
for the timing analysis. However, for the flux analysis, 16 phase bins were used across the 
pulse. For each phase bin, we maintained a spectral resolution of 128 bins over the PCA 
range. Given the broad morphologies of the average pulse profile, only one phase bin could 
be used as a background estimator. The pulse profiles were then phase aligned, so that 
the same off-pulse bin was used for background in every case. The remaining phase bins 
were summed, and their spectral bins regrouped using the FTOOL grppha, such that no bin 
had fewer than 20 counts after background subtraction. Energies below 2 keV and above 
10 keV were ignored, leaving 17-21 spectral channels for fitting, depending on the data set. 
The regrouped, phase-summed data sets, along with the background measurement, were 
used as input to the X-ray spectral fitting software package XSPEC 5 . Not all PCUs were on 
during our observations, but this was taken into account when producing response matrices. 
Response matrices were created using the FTOOLs xtef ilt and pcarsp. Because of the 
limited spectral resolution, fitting a two-component model was not practical, so we used a 
simple photoelectrically absorbed power-law, at first holding only Nh fixed (see Table 4). 
For all sources we found that the photon index V was constant within the uncertainties; we 
therefore held it fixed at its mean value. To extract a pulsed flux at each observing epoch 
we refitted each spectrum by varying only the normalization. Uncertainties were measured 
using the XSPEC command steppar. Our flux time series for 4U 0142+61, 1RXS 1708-4009 
and IE 2259.1+586 are displayed in Figures 9, 10 and 11 respectively. We do not find 
evidence for any large variability in the pulsed flux. Although the measured xt statistics for 
IE 2259.1+586 and 1RXS 1708-4009 are significantly larger than unity, given that the error 
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bars are statistical only, we cannot rule out unmodeled systematic errors. A long-term flux 
monitoring campaign with an imaging telescope is obviously desirable. 



4. Discussion 

4.1. Timing 

With our report of successful phase-coherent timing for 4U 0142+61, the detailed long- 
term timing properties of all known AXPs are becoming clear. As a population, they are 
generally spinning down steadily and with impressive stability. Indeed, for IE 2259.1+586, 
the stability is comparable to that of many radio pulsars: the Ag statistic (defined as the 
phase deviation due to v after 10 8 s; see Arzoumanian et al. 1994) is 0.43, consistent with 
what would be expected for a radio pulsar having a comparable spin-down rate, —0.78, given 
the large scatter in the radio pulsar values. For 4U 0142+61, fitting for v in either of the two 
best solutions yields \v\ — (5 ± 1) x 10~ 24 s~ 3 , implying A§ = 0.89, slightly larger than that 
expected for a radio pulsar with the same spin-down rate, —0.42, though still comparable 
given the scatter. For 1RXS 1708-4009, A 8 is significantly higher than would be predicted 
(2 versus 0.2), however this is consistent with glitch recovery. 

The stability of IE 2259.1+586 is particularly surprising given past reported large devi- 
ations from simple spin-down (e.g. Baykal & Swank 1996). This is in contrast to IE 1048.1- 
5937 which also showed apparent deviations from simple spin-down (see, e.g. Oosterbroek 
et al. 1998; Paul et al. 2000) but which we have verified in RXTE monitoring (Kaspi et al. 
2001). As argued by Baykal et al. (2000), it may be a coincidence that IE 2259.1+586 has 
shown great rotational stability during our observations only (but see also § 4.4). 

The current stability of IE 2259.1+586 allows for a test of the Melatos (1999) model of 
radiative precession. In that model, given well-defined but uncertain assumptions about the 
geometry and location of the magnetic field, a highly magnetized neutron star could deviate 
from sphericity and exhibit significant precession with period of a few years. This was 
suggested as an explanation for timing anomalies seen for IE 2259.1+586 and IE 1048.1- 
5937. The current data for IE 2259.1+586 rule out, with 99% confidence, an amplitude 
for precession 0.013 times that used by Melatos (1999) to explain the timing anomalies. 
Radiative precession was also ruled out by RXTE monitoring for IE 1048.1-5937 (Kaspi 
et al. 2001). Frequent large glitches in both these sources, as were suggested by Heyl & 
Hernquist (1999), also have not been confirmed. Also, neither 4U 0142+61 and IE 1841-045 
(Gotthelf et al. 2001) show any evidence for either radiative precession or frequent large 
glitches. 
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Overall, however, a striking property of the timing of the AXPs as a whole, apart from 
the stability, is the diversity of behavior. IE 2259.1+586 and 4U 0142+61 have shown great 
stability; IE 1841-045 shows stability but with considerably more red noise (Gotthelf et al. 
2001); 1RXS 1708-4009 has exhibited a glitch with recovery similar to what is observed for 
radio pulsars; and IE 1048.1-5937 is less stable than even some known accretors (Kaspi et al. 
2001). 

Interestingly, there is a correlation between timing stability and v in our data. The 
sources with by far the smallest z>'s, IE 2259.1+586 and 4U 0142+61, are the most stable 
(at least during our observations), while that with the largest, IE 1048.1-5937, is the least 
so. The SGRs, being even less stable (Woods et al. 2000; Woods et al. 2001), have even 
larger z>'s, in agreement with this trend. IE 1841-045 and 1RXS 1708-4009 have timing 
stabilities and v intermediate between the extremes. One therefore might expect IE 1841- 
045 to exhibit a glitch in the near future. If correct, the trend suggests a continuum of 
timing properties between the AXP and SGR populations, lending additional support to the 
connection between them. Furthermore, it has long been recognized that radio pulsar timing 
stability is also correlated with spin-down rate (e.g. Cordes & Helfand 1980; Arzoumanian 
et al. 1994). This suggests, in addition, a connection between all three populations, and a 
continuum of timing behavior that depends on stellar magnetic field. If correct, this provides 
strong support for the magnetar model. 

Marsden & White (2001) noted a correlation between spin-down rate and spectral prop- 
erties in AXPs. An additional correlation with timing stability could strengthen the argu- 
ment that the spin properties, emission mechanisms, and internal structure and dynamics 
are all related in these objects. 



4.2. Pulse Morphology 

The presence of substantial harmonic content in AXP pulse profiles (see Fig. 8) is at 
odds with models in which the X-ray emission is thermal from the surface of a hot neutron 
star having conventional magnetic field. As shown by Dedeo et al. (2001), the expected 
range of ratios of second to first Fourier amplitudes, A2/A1, can be predicted for a neutron 
star emitting thermal X-rays, as a function of total pulse fraction, for different radially 
peaked beaming functions, surface hot-spot geometries, and neutron star compactnesses. 
Low harmonic content is expected primarily because of gravitational light bending, unless 
there is significant beaming. Given the pulsed fraction of IE 2259.1+586 (Ozel et al. 2001, 
and references therein), Dedeo et al. (2001) show that the ratio of A2/A1 < 1 for all possible 
geometries and compactnesses, while Figure 8C shows A 2 /Ai = 3.97 ± 0.01 in the soft 2- 
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4 keV band. The presence of a significant third harmonic in 1RXS f 708-4009 (Fig. 8B) and 
even a significant sixth harmonic in IE 2259. 1 +586 (Fig. 8C) is similarly problematic. On 
the other hand, Ozel (2001) suggests that, in the magnetar model, when properly accounting 
for the propagation of X-rays through the neutron star atmosphere, the pulse morphology 
can be surprisingly non-sinusoidal. In particular, in 10 14 — 10 15 G magnetic fields, the 
emission from a surface patch can be highly and non-radially beamed. Furthermore, if 
the surface distribution itself is non-isotropic, AXP profiles may be reproduceable. Indeed, 
Ozel, Psaltis & Kaspi (2001) show that AXP pulsed fractions and luminosities can only be 
reproduced, assuming thermal emission from the surface of a magnetar propagating through 
a hydrogen atmosphere, with a single hot spot. In addition, Thompson, Lyutikov & Kulkarni 
(2001) argue that AXP pulse profiles can be accounted for even if the surface emission were 
isotropic, because of resonant scattering off charged particles in the magnetosphere, which 
has anisotropic optical depth. A detailed comparison of model predictions with observed 
AXP pulse profiles may prove interesting. 



4.3. Pulse Morphology Energy Dependence 

The energy dependence of the pulse profiles of AXPs can in principle be a strong con- 
straint on the emission model. We have studied this by comparing the morphologies of 
pulse profiles in two energy bands. The AXPs display a variety of behaviors in this regard. 
For example, while IE 1048.1-5937 and IE 1841-045 show essentially no change in pulse 
morphology with energy (Kaspi et al. 2001; Gotthelf et al. 2001), 1RXS 1708-4009 shows 
tremendous variation (Fig. 7B). In no case (except possibly 1RXS 1708-4009) are the pro- 
files in the two energy bands very different, problematic for interpreting AXP two-component 
spectra as being a result of two independent emission mechanisms (e.g. Perna et al. 2001). 
Figure 7A shows a possible energy-dependent phase lag, in which the low energy pulse lags 
the high energy pulse, however, this is not observed in any other AXP. 

There is a possible rough trend in which sources with the least harmonic content 
(IE 1048.1-5937, IE 1841-045) show the least energy dependence. However the lack of strong 
pulse morphology evolution with energy in IE 2259.1+586 compared with 1RXS 1708-4009 
does not obviously support this. Another possibly interesting trend is that for the sources 
in which the ratio of the Fourier amplitudes of the second and first harmonics is large 
(IE 2259.1+586, 4U 0142+61), the ratio is largest at low energies. This manifests itself in 
the profile plots as the ratio of peak heights being closer to unity in the hard band than in 
the soft, which is true of 1RXS 1708-4009 as well. 

Whether and how these observations can constrain models remains to be seen. Ozel 
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(2001) shows that, when photon propagation through a magnetar atmosphere is properly 
modeled, the angular dependence of the emission has two components: a narrow pencil beam 
at small angles with respect to the surface normal, and a broad fan beam at large angles. She 
shows that the importance and opening angle of the pencil beam increases with increasing 
photon energy. This may be able explain the observed change in the relative peak amplitudes 
as a function of photon energy, especially in 1RXS 1708-4009, if part of the pulse is from a 
fan beam and part from a pencil beam. However detailed pulse profile modeling is required. 



4.4. Pulsed Flux Time Series 

The pulsed fluxes of all known AXPs are roughly constant, with RMS fluctuations of 
at most ~20-30% of the mean (la) (see Table 4, Kaspi et al. 2001, Gotthelf et al. 2001). 
Past studies of the flux stability of two AXPs, IE 2259.1+586 and IE 1048.1-5937, have 
suggested that their phase-averaged fluxes are highly variable, with fluctuations possibly 
as large as a factor of ten (Baykal & Swank 1996; Oosterbroek et al. 1998). Large flux 
variations, particularly on time scales under a year and in the absence of torque variations, 
are unexpected in the magnetar model. This is because in this model, X-rays come from the 
hot surface where the heat is maintained by the decay of the interior magnetic field. As the 
thermal conduction time from core to surface is roughly a year (Van Riper et al. 1991), flux 
variations on shorter time scales are hard to explain in the absence of burst activity as is 
seen in the SGRs (e.g. Woods et al. 2001). 

With the caveat that we measure pulsed flux while other studies measured total flux, it 
is surprising that we find such flux stability in our observations given past claims. Baykal 
et al. (2000) argued that the flux stability of IE 2259.1+586 is a result of its currently 
being in a quiescent accreting state that is also characterized by timing stability (see also 
§ 4.1). Though we cannot rule this out, we point out that (i) the same coincidence would 
have to be true of IE 1048.1-5937 for flux stability, (ii) the pulsed fluxes of 4U 0142+61, 
1RXS 1708-4009 and IE 1841-045 are also stable, and (iii) the pulsed flux of IE 1048.1-5937 
is stable even though it shows unstable timing behavior, which argues that the flux stability is 
independent of spin-down behavior. We note that our flux measurements, in contrast to past 
studies, have been made using a single instrument, single bandpass, and analysis method, 
which eliminates the difficulty in comparing measurements made with different instruments, 
different spectral ranges, and different analysis methods. 

Vasisht et al. (2000) have suggested that the variable X-ray pulsar AX J1845. 0-0258 
which is at the center of the supernova remnant G29. 6+0.1 may also be an AXP. This source 
has shown variations in its X-ray luminosity by a factor of ~10 between observations taken 
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6 yr apart. Given that the pulsed fluxes of all known AXPs are stable, the interpretation of 
AX J1845. 0-0258 as a bona fide AXP is therefore questionable. 



5. Conclusions 

We have presented a variety of observational results for three anomalous X-ray pulsars, 
4U 0142+61, 1RXS 1708-4009 and IE 2259.1+586, obtained using regular monitoring obser- 
vations and archival observations from the PCA aboard RXTE. Our results, combined with 
those for IE 1048.1-5937 and IE 1841-045 (Kaspi et al. 2001; Gotthelf et al. 2001), provide 
significantly improved descriptions of AXP X-ray properties that can be used to test models 
for the nature of these sources. 

Specifically, with our successful phase-coherent timing of 4U 0142+61, the detailed tim- 
ing properties of all known AXPs are now established and can be compared. We find a wide 
variety of timing behaviors, ranging from high stability (in IE 2259.1+586 and 4U 0142+61), 
to instabilities so severe that phase-coherent timing is not possible (in IE 1048.1-5937). We 
note that timing stability appears correlated with decreasing v. If this trend is real, it sug- 
gests a continuum of timing properties between the AXP and SGR populations, lending 
additional support to the connection between them. We note that a correlation between 
timing instability and spin-down rate has long been recognized among radio pulsars (e.g. 
Cordes & Helfand 1980; Arzoumanian et al. 1994), suggesting a connection between all three 
populations. 

We have also used the RXTE data to investigate a variety of other AXP properties. 
Motivated by a report of a significant change in pulse morphology in IE 2259.1+586 in 1990 
(Iwasawa et al. 1992), we searched for changes in the pulse morphology at each epoch for all 
targets. We find no significant pulse morphology variations. Typically, we rule out variations 
in features having amplitude >20% of the peak amplitude at the la level, although the limit 
depends on source and integration time. 

We have presented high signal-to-noise average pulse profiles for each AXP, and con- 
sidered them as a function of energy. We show that, as in the timing properties, there is a 
variety of different behaviors for the energy dependence. Possible trends include a greater 
energy dependence for pulses profiles having greater harmonic content, and, for the latter 
sources, relative peak amplitudes tending closer to unity as photon energy increases. De- 
tailed modeling of AXP pulse profiles, in addition to their pulsed fractions and luminosities 
(Ozel et al. 2001) can test models in which the emission is from the surface of a cooling 
magnetar (Ozel 2001; Ho & Lai 2001; Zane et al. 2001). 
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Finally, we use the monitoring and archival data to obtain pulsed flux time series for each 
source. We have found no large changes in pulsed flux for any source, and have set la upper 
limits on variations ~20-30% (depending on the source). This is surprising given previous 
reports of large (factor of 5-10) total flux variations in IE 2259.1+586 and IE 1048.1-5937 
(Baykal & Swank 1996; Oosterbroek et al. 1998). Assuming a constant pulsed fraction, this 
suggests that more than one of the AXPs happen to be much more quiescent during the 
RXTE monitoring than in the past. 
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Table 1. Summary of RXTE Observations 



Observing 


Period 


Nominal Nominal 


Number 


Total 






Exposure Separation 


of Obs. 


Exposure 






(ksec) (weeks) 




(ksec) 




4U 0142+61 


Mar 1996 






1 


7.0 


Mar 1996 






1 


43.5 


Mar 1998 






1 


20.3 


Nov 1996 - 


- Dec 1997 


1.1 4-5 


14 


14.9 


Mar 2000 


- Feb 2001 


3.3 4-5 


14 


46.0 


Mar 2001 


- Jun 2001 


6.6 4-5 


3 


19.9 


1RXS 1708-4009 


Jan 1998 - 


Jan 1999 


6.2 4-5 


12 


74.0 


Feb 1999 - 


Jun 2001 


3.1 4-5 


33 


103.6 


May 2001 




2.0 1/7 


5 


10.1 


IE 2259.1+586 


Sep 1996 






1 


86.9 


Nov 1996 


- Dec 1997 


1.0 4-5 


14 


14.0 


Feb 1997 - 


Mar 1997 


24.8 1 


5 


124.1 


Aug 1998 


- Sep 1998 


14.0 1 


8 


111.7 


Jan 1999 - 


Mar 2001 


3.2 4-5 


19 


61.3 


Mar 2000 


- May 2001 


7.0 4-5 


11 


77.3 
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Table 2. Spin Parameters for 4U 0142+61. f 



First Span 



Second Span 



A 



B 



MJD Range 
No. TOA 
v (Hz) 

v (10~ 14 Hz/s) 
Epoch (MJD) 
RMS residual 



50411-50893 51610-52028 50411-52028 50411-52028 

15 17 32 32 

0.115096877(10) 0.1150969336(5) 0.1150969299(10) 0.1150969209(11) 

-2.687(5) -2.649(11) -2.5980(23) -2.5969(24) 

51704.000 51704.000 51704.000 51704.000 

0.027 0.011 0.074 0.077 



a Ephemerides A and B are both viable given the data set, which includes a 2-yr gap at 
the center of the 4.4-yr span. Ephemeris B has one additional pulse in the gap relative to 
Ephemeris A. See § 3.1. 
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Table 3. Spin Parameters for 1RXS 1708-4009 and 
IE 2259.1+586. 





1RXS 1708-4009 


IE 2259.1+586 


MJD Range 


51472-51995 


50356-52016 


No. TOA 


19 


67 


v (Hz) 


0.090917063(5) 


0.14328806234(8) 


v (10~ 14 Hz/s) 


-16.07(2) 


-0.99434(16) 


v (10- 23 Hz/s 2 ) 


5.2(6) 


0.228(14) 


Epoch (MJD) 


51215.9308 


51995.5827 


RMS residual (periods) 


0.0075 


0.0099 
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Table 4. AXP Spectral Parameters and 2-10 keV Pulsed Flux Variability Upper 

Limits 





Spectral 


Parameters 


Flux Time Series 


Parameters 


Source 


\ogN H a 
(cm~ 2 ) 


Refs. b T c 


xl d Degrees of 
Freedom 


RMS flux c 
mean flux 


4U 0142+61 


22.18 


1 3.50 


0.63 36 


0.26 


1RXS 1708-4009 


22.25 


2 2.63 


1.45 52 


0.17 


IE 2259.1+586 


22.08 


3 4.02 


1.38 59 


0.28 



a For each AXP, Nh was held fixed at these values. 

b (l) Mereghetti & Stella (1995); (2) Sugizaki et al. (1997); (3) Baykal & Swank 
(1996). 

c Mean power-law index as measured for the pulsed portion of the spectrum. 

d xl statistic of the pulsed flux time series displayed in Figures 9, 10 and 11. 

e RMS of the pulsed flux time series displayed in Figures 9, 10 and 11 in units of 
the mean pulsed flux. 
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Fig. 1. — Average pulse profiles of the AXPs. Two cyles are plotted for clarity. 
(A) 4U 0142+61 (Total exposure time: 149 ks, energy range: 2.5 - 9.0 keV); (B) 1RXS 1708- 
4009 (185 ks, 2.5 - 9.0 keV); (C) IE 2259.1+586 (469 ks, 2.2 - 5.5 keV). 
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Fig. 2. — Harmonic content of the average pulse profiles in Figure 1. (A) 4U 0f42+61 
(Energy range: 2.5 - 9.0 keV) ; (B) fRXS f 708-4009 (2.5 - 9.0 keV); (C) IE 2259. f +586 
(2.2 - 5.5 keV). The ratio of the power in the n harmonic to the total power in all harmonics 
is plotted versus n. 
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Fig. 3. — Arrival time residuals for 4U 0142+61 with Ephemeris A (Table 2), that includes 
v and z>, subtracted. Residuals relative to Ephemeris B are qualitatively similar. 
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Fig. 4. — Pulse frequency history for 4U 0142+61. Top panel: the solid, thick line is the 
ephemeris detemined from our phase- coherent timing with RXTE. Individual data points 
are from a variety of X-ray missions (see Israel et al. 1999b for references) and are shown 
with 90% confidence level error bars. Note that on this scale, Ephemeris A and B (see § 3.1, 
Table 2) are indistinguishable. The dotted line is the extrapolation of our fit backward. 
Bottom panel: the same data but with the linear trend and offset removed. In both plots, 
the uncertainties (including that due to the ambiguity between A and B) are smaller than 
the width of the lines. 
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Fig. 5. — Arrival time residuals for 1RXS 1708-4009 with u, v and v subtracted. 
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Fig. 6. — Arrival time residuals for IE 2259.1+586 with v, v and i) subtracted. 
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Fig. 7.- Average pulse profile of (A) 4U 0142+61, (B) 1RXS 1708-4009 and 
(C) IE 2259.1+586 in two energy bands as observed by RXTE. Two cycles are plotted 
for clarity. Note that the scaling was chosen to minimize the x 2 of the difference between the 
two profiles. Thus the only information that these plots convey is the relative amplitudes of 
the features of the profile. 
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Fig. 8. — Harmonic content of the average pulse profiles in Figure 7 in three energy bands. 

(A) 4U 0142+61; (B) 1RXS 1708-4009; (C) IE 2259.1+586. The ratio of the power in the 

th 
n harmonic to the total power in all harmonics is plotted versus n. 
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Fig. 9. — Pulsed flux time series for RXTE observations of 4U 0142+61. Error bars represent 
la confidence intervals. See §3.4 for details of the analysis procedure. 
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Fig. 10. — Pulsed flux time series for RXTE observations of 1RXS 1708-4009. Error bars 
represent la confidence intervals. See §3.4 for details of the analysis procedure. 
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Fig. 11. — Pulsed flux time series for RXTE observations of IE 2259.1+586. Error bars 
represent lcr confidence intervals. See §3.4 for details of the analysis procedure. 



